Nonlinear dynamics and chaos in parametric sound generation 



Victor J. Sanchez-Morcillo, Victor Espinosa, Javier Redondo, and Jesus Alba 
Departament de Fisica Aplicada, Universitat Politecnica de Valencia, 
Crta. Natzaret-Oliva s/n, 4.6730 Grau de Gandia, Spain 

Short title: Chaos in acoustic interferometers and V. J. Sanchez-Morcillo, V. Espinosa, J. Redondo, and J. Alba 
Departament de Fisica Aplicada, Universitat Politecnica de Valencia, 
Crta. Natzaret-Oliva s/n, 46730 Grau de Gandia, Spain 

Abstract 

A theoretical analysis of the subharmonic generation process in an acoustical resonator (interferometer) with plane walls is 
performed. It is shown that, when both the pumping wave and the generated subharmonic are detuned with respect to the 
resonator modes, the fields can display complex temporal behaviour such as self-pulsing and chaos. A discussion about the 
acoustical parameters required for the experimental observation of the phenomenon is given. 



PACS numbers: 43.25.Ts, 43.25.Rq 
I. INTRODUCTION. 

Nonlinear dynamical phenomena have been reported 
in several acoustical systems. Different complex scenar- 
ios of temporal evolution, including self-pulsing, period 
doubling and deterministic chaos, are typically observed 
in a large variety of driven oscillators, such as nonlinear 
resonators, musical instruments or ultrasonic cavitation 
(for a review, see Ref. 1). In the particular resonator 
configuration, self-action of sound propagating in a vis- 
cid fliud, where the thermal mechanism of nonlinearity is 
dominant, has been studied theoretical and experimen- 
tally, and the onset of bistable, self-pulsed and chaotic 
regimes has been shown to occur with precedence to an- 
other acoustical nonlinear effects related to the elastic 
(quadratic) nonlinearity of the medium. These non- 
thermal effects, involving the appearance of new frequen- 
cies not present in the driving source, provoke the exci- 
tation of higher harmonics of the driving source (which 
in the absence of dispersion leads to waveform distortion 
and ultimately to the development of shock waves) or the 
frequency division of the input wave (subharmonics), an 
effect so called parametric sound generation. The latter 
phenomenon consist in the resonant interaction of a triad 
of waves with frequencies u)q,u>i and u>2, which are re- 
lated by the energy conservation condition uq = lo\ + u>2. 
The process is initiated by an input pumping wave of 
frequency loq which, due to the propagation in the non- 
linear medium, generates a pair of waves with frequen- 
cies oj\ and u>2- When the wave interaction occurs in a 
resonator, a threshold value for the input amplitude is 
required. This process has been described before by sev- 
eral authors under different conditions, either theoretical 
and experimentally. 0, 0, IE @ 

It is important to note that the presence of higher har- 
monics can be avoided by different means. One method 
considers that the interaction takes place in a dispersive 
system. Dispersion allows that only few modes, those sat- 
isfying given synchronism conditions, participate in the 
interaction process. In finite geometries, such as waveg- 



uides or resonators, the dispersion is introduced by 
the lateral boundaries. Different resonance modes, prop- 
agating at different angles, propagate with different ef- 
fective phase velocities. However, in homogeneous un- 
bounded systems dispersion is usually not present. Dif- 
ferent dispersion mechanisms have been proposed in non- 
linear acoustics, such as propagation in bubbly or layered 
(periodic) media, which makes sound velocity prop- 
agation to be wavelength dependent. Other proposed 
mechanisms are, for example, media with selective ab- 
sorption, in which selected spectral components experi- 
ence strong losses and may be removed from the wave 
field. , or resonators where the end walls present some 
frequency-dependent complex impedance. In this case, 
the resonance modes of the resonator are not integrally 
related, and by proper adjustment of the resonator pa- 
rameters one can get that only few modes, those lying 
close to a cavity resonance, reach a signifficant ampli- 
tude. Under this conditions, a spectral approach to the 
problem results appropriate, as has been proposed and 
experimentally confirmed in Ref. 6. 

Therefore the selective effect of the resonator allows to 
reduce the study of parametric sound generation to the 
interaction of three field modes, correspondig to the driv- 
ing (fundamental) and subharmonic frequencies, and to 
describe this interaction through a small set of nonlinear 
coupled differential equations. In the present work we 
describe the particular degenerate case of subharmonic 
generation, where wi = 102 and consequently luq = 2uj\, 
being ujq the fundamental and lo\ generated subhar- 
monic, both quasi-resonant with a corresponding reso- 
nance mode. This degenerate case has been the matter of 
previous experimental studies, 0, Q and the observation 
of dynamical behaviour has been reported. 0, ^| How- 
ever, despite the experimental observation of complex 
temporal dynamics in this system, a theoretical frame- 
work supporting these effects and taking into account the 
peculiarities of the acoustic systems is absent. The main 
purpose of this work is to establish the necessary con- 
ditions for the development of nonlinear dynamics in an 



acoustic resonator where parametric generation of sound 
takes place. 



II. MODEL EQUATIONS. 

We consider a resonator with plane walls with finite 
width (acoustic interferometer) driven at a frequency u>q, 
and that the wave interaction is only effective among two 
resonant frequencies: the pump loq and the subharmonic 
u>i, for which the relation ujq = 2u>\ hold. The field inside 
the resonator can then be expanded as 



'(r,*) = £) P ^ r '*)' 



(1) 



where Pj (r, i) are the wave components related with 
the frequency uij. Restricting the analysis to a one- 
dimensional resonator, where quasi-plane waves propa- 
gate on the z axis, and taking into account that, due 
to reflections in the walls, there exist waves propagating 
simultaneously in opposite directions, the wave compo- 
nents take the form 



Pj (z,t) = pj (t) cos (kjz) 



e 3 + c.c. 



(2) 



where fc| a cavity eigenmode. 

Under the assumption of plane waves inside a res- 
onator with highly reflecting ends (weak cavity losses), 
the model for the time evolution of degenerate paramet- 
ric sound generation reads [l2| 
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E - 7o (l + iA Q ) P() - i/3pl, 
-7 1 (l+iA x )p 1 -i/3p*p , 
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where the variables pi (i = 0, 1) are the complex slow- 
varying amplitudes of the fundamental and subharmonic 
fields, with angular frequencies Ui. The rest of parameters 
and their units are: the pump E (Pa/s), proportional to 
the injected amplitude pi n {Pa) of the fundamental mode 
and acting as a control parameter; the loss coefficient 7, 
(s _1 ); the nonlinear coupling parameter /3 (s _1 _Pa _1 ), 
and the adimensional detuning parameter A,, propor- 
tional to the frequency mistuning between the field fre- 
quency LOi and the closest resonator eigenfrequencies w". 
All these parameters are defined as 

E = -^j-\/%Pin, 



cT 



A, = 
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71 



(4) 
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being c the propagation velocities, % and IZi, the 
transmission and reflection coefficients at resonator ends 



(IZi » %, has been assumed), L the length of the res- 
onator, £ the nonlinearity parameter, p the medium den- 
sity and Ki — 5kf/2 account for losses related to diffusiv- 
ity of sound in the medium 8. A detailed derivation of 
this model, generalized to include the field dependence 
on the transversal coordinates (not considered here for 
simplicity) can be found in Ref. 15. The same equations 
have been obtained in Ref. 3 for waveguide resonators 
(althought with different coupling coefficients, dependent 
on the structure of the interacting transverse modes), or 
for one-dimensional acoustic interferometers in Ref. 6, in 
a form identical to Eqs. ©■ 

The model (01 has been also derived in the context of 
nonlinear optics. |13( where the existence of complex tem- 
poral dynamics has been reported. In the next section 
we review the main results concerning the solutions and 
its stability properties for our particular acoustic system. 



III. STATIONARY SOLUTIONS. 

Two stationary solutions of Q are obtained when the 
temporal derivatives vanish. || The simplest case corre- 
sponds to the trivial solution, 



E 



Po 



7o(l + *A ) 



, Pi = 0, 



(5) 



characterized by a null value of the subharmonic field 
inside the resonator. Increasing the pump value E this 
solution bifurcates into a new solution, in which also the 
subharmonic field has a nonzero amplitude \p 1 \ , given by 



IPil 2 = i(7o7i(A A 1 -l)± 



(6) 



V //3 2 |i?| 2 -7g7 1 (Ao + AO 2 ), 



while the stationary fundamental amplitude \p \ takes 
the value 



\Po\ = j\/l + A?, 



(7) 



The emergence of this finite amplitude solution corre- 
sponds to the process of subharmonic generation. Note 
that the fundamental amplitude above the threshold is 
independent of the value of the injected pump, which 
means that all the energy is transferred to the subhar- 
monic wave. 

The transition from solutions 10 to © occurs at a 
critical pump amplitude 



\E\ 



(3 



^/(1 + AgHl + A?), 



(8) 



beyond which the trivial solution is unstable. These 
results have been confirmed experimentally in Ref. 3. 
The character of the bifurcation depends on the detun- 
ing values. As demonstrated in Ref. 3, and also in the 



2 



optical context, 01 the bifurcation is supercritical when 
A0A1 < 1, and subcritical when A0A1 > 1. In the latter 
case, both trivial and finite amplitude solutions can co- 
exist for given sets of the parameters, which results in a 
regime of bistability between different solutions. Figure 1 
illustrates the form of the solutions in these two regimes. 
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FIG. 1: Bifurcation diagram of the subharmonic field in the 
monostable (a) and bistable (b) cases. Subharmonic detuning 
values are Ai = 2.4 (a), and Ai = —2.4 (b). The rest of pa- 
rameters are Ao = 1, j3 = 0.002 and 7o = 7i = 500 (see text 
below). Pressure units are given in Pa. Unstable branches are 
depicted in dashed lines, and the black square correspond to 
a Hopf bifurcation point. 



IV. LINEAR STABILITY ANALYSIS. 

In order to demonstrate the possibility of temporal dy- 
namics in this system, the stability of the previous sta- 
tionary solutions has been investigated by substituting 
in Eqs. © and their complex conjugate a perturbed so- 
lution in the form 



Pi (t) = Pi + Spi(t) =Pi + 



-At 



(9) 



where p i are the particular stationary values. By lin- 
earizing the resulting equations around the small per- 
turbations 5pi, one obtains A as the eigenvalues of the 
stability matrix that relates the vector of the perturba- 
tions (Spi(t), Sp*(t)) with their temporal derivatives. The 
instability of the solution is determined by the existence 
of positive real parts of the roots A of the characteristic 
polynomial, 



where, for the finite amplitude solution Eq. JBJ) and J7J, 

the coefficients take the form 



ci = 2( 7o + 7l ), 

c 2 =4/? 2 ||5 1 | 2 +7o + 47o7i+7oA§, 

c 3 = 4/3 2 1^1 2 ( 7o + 7i) + 2 7 o7i (l + A 2 ) 

C4 = 4/3 2 Ip^ 2 |p 1 | 2 -7 07l (AoAi-i; 



(11) 



Instead of solving the quartic polynomial Eq. (|10|) to 
find the eigenvalues, we apply the Hurwitz criterion (see 
e.g. Ref. 17), which states that the solution becomes 
unstable when at least one of the following conditions 
are violated: 



ci > 0, c 2 ;> 0, c 3 > 0, c 4 > 

cic 2 - c 3 > 

C3 (cic 2 - c 3 ) - c 2 c 4 Js 0. 



(12a) 
(12b) 
(12c) 



Condition (|12bjl is always fulfilled. In conditions (|12afl . 
only C4 can become negative, corresponding with the 
turning point of the stationary solution of the A Ai > 
1 case (see Fig. lb). This means that the lower branch 
of the stationary solution, represented by the dashed line 
in Fig. 1 is always unstable. 

Finally, the equality in condition 112cfl denotes the ex- 
istence of a pair of complex conjugate eigenvalues with 
vanishing real part (A = Hlj), which are indicative of a 
Hopf bifurcation which gives rise to oscillating solutions, 
with angular frequency u>. Substituting (|ll|l in (|12cfl the 
instability condition reads 

7^1(1 + A 2 ) ((7 + 2 7l ) 2 +7o : A 2 ) + (13) 

2 (7o + 7i) 2 (7o (1 + Ao) + 2 7l (1 + AoAr)) (3 2 \p, | 2 ^ 0. 

which is fullfilled when the subharmonic amplitude 
reaches the value 



I- I 2 

\Pi\hb = 



-7§7i (1 + Aj) ((7 + 2 7l ) 2 +7gA 2 ) 

W 2 (7o + 7i) 2 (7o (1 + Al) + 2 7l (1 + AqAO) ' 

(14) 



The positiveness of the threshold value l|14|) requires an 
additional condition, namely that the product of the de- 
tunings satisfies 



AoAx < --L( 7o + 2 7l + 7o A 2 ) 
z 7i 



(15) 



A + ciA J + c 2 A 2 + c 3 A + c 4 = 



(10) 



The latter condition excludes bistability (which re- 
quires A0A1 > 1). Consequently, condition 115|) implies 
that the bifurcation that leads to subharmonic genera- 
tion is always supercritical, and then only single-valued 
solutions can undergo a Hopf bifurcation (see Fig. 1). 
The angular frequency of the oscillations at the Hopf in- 
stability point can be found by substituting A = iuj in 
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(|13fl . and reads 

,, - htfZ 2 i ^i( 1 + a q) na\ 

V 7o + 7i 

where is the subharmonic amplitude at the Hopf 

bifurcation point given by Eq. i|14[l . 

V. ACOUSTICAL ESTIMATES. 

For a better understanding of the predictions of the 
model and to assure their experimental validity, a special 
care must be taken when choosing the acoustical parame- 
ters involved in the theory. We will first refer such quan- 
tities to the experimental investigation of subharmonic 
generation in an acoustic interferometer described in de- 
tail by Yen [|| . It consists in a three sectioned resonator 
(quartz crystal— water— quartz crystal) of variable length, 
operated at 26° C. Our aim is to evaluate the threshold 
for subharmonic generation Eq. @ and the possibility of 
its destabilization at the value given by l|14|). From Ref. 
6 we obtain realistic acoustical values for the field decay 
rates r y i , the detunings A, and the nonlinearity parame- 
ter (3 involved in the model. These parameters depend, 
as defined in J3J|, on sound frequencies and velocities, 
medium nonlinearity s, cavity length and reflectivity of 
boundaries, and the eigcnfrcqucncics of the resonator. In 
Ref. 6 the resonance modes of the one-dimensional res- 
onator were calculated by using the formula 

R(t&nk' w D + tank' w H) +tanfc'5- (17) 
R 2 tan k' w D tan k' w H tan k'S = 

where R is the ratio of wall impedance to medium 
impedance, k' and k' w are the wavenumbers in the fluid 
and the walls respectively, and D, H and S are the thick- 
nesess of the two end walls and the length of the ac- 
tual resonator, respectively. From l|17|l results an non- 
equidistant spectrum, which is required to allow the 
fulfillment of the condition l|15|) . In the case of an 
equidistant spectrum (corresponding to a perfect res- 
onator without losses), the detunings are obeyed to have 
the same sign. 

The values of these resonance modes, as well as the 
quality factor Qi associated with them, were also exper- 
imentally measured in Ref. 6. The pumping frequency 
was about 1.6 MHz, for which a subharmonic emission 
of 0.8 MHz is expected. The closest resonator eigcnfrc- 
quencies to such fundamental and subharmonic emissions 
were / C = 1.5853 MHz and /f = 0.8088 MHz, with asso- 
ciated quality factors of 4700 and 2600, respectively. The 
quality factor allows to estimate directly the decay rates 
7j using the relation Qi = uJi/2^. h and the corresponding- 
detuning is calculated then from . The nonlinear cou- 
pling parameter (3 defined in Eq. is determined by 
the subharmonic frequency u>±, the medium nonlinearity 
(e = 3.5 for water), and the sound velocity c = 1500 m/s. 



Note that the possible values of the detunings Ao and Ai 
are restricted by the non-equidistant eigenfrequencies of 
the resonator and the parametric character of the pro- 
cess, since ujq = 2u>i. 

The threshold of subharmonic generation was also 
measured by Yen,|(| and given in terms of the funda- 
mental pressure amplitude inside the resonator at the 
bifurcation point given by Eq. J7J. The model predicts 
its dependence on the subharmonic field detuning and 
on the subharmonic frequency through the coupling pa- 
rameter (3. The lowest theoretical value of the threshold 
corresponds to a perfect tuning of the subharmonic field 
with a cavity eigenmode, that is for a fixed fi — 0.8088 
MHz and the corresponding fundamental fo = 1.6176 
MHz, results Ai = 0, and Ao = —187.7, and nonlinear 
coupling parameter (3 — 0.002 s _1 Pa~ 1 . The predicted 
fundamental pressure amplitude at threshold is therefore 
about [po| = 0.49 MPa or 4.9 bar. The lowest threshold 
reported in Ref. 6 is about 8.8 bar, measured for a fre- 
quency close to the cavity mode /q = 1.5853 MHz. This 
difference can be due to the finite bandwidth of the excit- 
ing transducer in opposition to the monochromatic pump 
considered here, since the introduction of small detunings 
on the expected subharmonic field causes a sensitive and 
rapid increase of the subharmonic emission threshold. 

After the switching-on of the subharmonic field the 
system supports constant pressure amplitudes of both 
fundamental and subharmonic fields, corresponding to 
the solutions ijfjj and J7J. In order to observe self- 
oscillations of the amplitudes, the Hopf bifurcation con- 
dition Q15f) must be fulfilled, e.g. by increasing the de- 
tuning of the subharmonic field. One would get this 
by tuning the frequency of the subharmonic up to the 
value /i = 0.8072 MHz, for which Ai = 52.4551 and 
Ao = —94.7686. Unfortunately, when trying to achieve 
the detunings condition for instabilities, the subharmonic 
emission threshold grows up to the unrealistic value 
|jPol = 26 MPa. For the same parameters, the Hopf bi- 
furcation threshold is obtained for a subharmonic ampli- 
tude of \Pihb\ = 1-14 GPa. Clearly, for this particular 
resonator mode distribution the operation regime should 
be always in the steady state with constant pressure am- 
plitudes. This could be the reason why dynamical states 
have not been described neither in Ref. 6 nor in other ex- 
perimental works in which the conditions were probably 
similar. 

Nevertheless the destabilisation of the subharmonic 
regime would explain certain phenomena previously de- 
scribed in the literature, such as self-modulation of the 
fields or possible chaotic behaviour |n|]. We esti- 
mate next the necessary conditions for the occurrence 
of these phenomena at reasonable values. The resonator 
described in Ref. 11 corresponds to that described above; 
two air backed quartz transducers placed in a Fabry- 
Perot configuration with water, with a regulable separa- 
tion between 5 and 10 cm. The driven transducer had a 
resonant frequency /o = 1 MHz and the other had a lower 
resonance frequency. We take this frequency as the fun- 
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damental for our estimations and evaluate the remaining 
parameters of the desired resonator. From inspection of 
Eq. (jHl it follows that the subharmonic emission thresh- 
old can be lowered by improving the quality of the res- 
onator for both modes. However, the behaviour of the 
second threshold it is not so evident: Fig. 2 shows the 
dependence of the pressure amplitude of the subharmonic 
field at the Hopf bifurcation on the ratio of the quality 
factors affecting to both fundamental and subharmonic 
fields. It evidences that instabilities are reached earlier 
when the fundamental field has significant lower losses 
than the corresponding subharmonic. According to this, 
we set Qo = 6000 and Qi = 3000, which are realistic 
values if compared to those measured in Ref. 6. 





Qo/Qi 



FIG. 3: Detuning combinations required to fulfil the Hopf 
bifurcation condition. 



The linear stability analysis can not predict the be- 
haviour of the fields beyond the bifurcation points, and 
also does not give information about the subcritical or 
supercritical character of the bifurcation. In the next sec- 
tion we perform the numerical integration of the model 
(0 using the parameters estimated above. 



VI. NUMERICAL SIMULATIONS. 



FIG. 2: Dependence of the pressure amplitude of the subhar- 
monic field at the Hopf bifurcation (\Pi\ HB ) on the ratio of the 
quality factors of fundamental and subharmonic fields. The 
parameters are A = 1, Ai = -2.4, j3 = 0.002 and Qi = 3000. 
The value of Qo is kept free. 

The resonator must meet another important restriction 
concerning its frequency distribution, in order to satisfy 
Eq. (|15fl . In Fig. 3 are shown the combination of detun- 
ings that accomplish this condition for the given quality 
factors. The shadowed area represent detuning values 
which are compatible with the Hopf bifurcation condi- 
tion (the inverted symmetrical figure, with opposite signs 
of both detunings, is also valid). One must also consider 
the lowest possible detunings in order to achieve the min- 
imal value for the pressure amplitudes in both thresholds. 
Note that, whilst Ao can be arbitrarily small, the abso- 
lute value of Ai must be greater than 2. A combination 
of values that reduce the thresholds to common exper- 
imental levels is (Ao,Ai) = (1.0,-2.4), corresponding 
to an intracavity pressure amplitude of the fundamental 
mode at the first threshold of |JJ | = 1.11 MPa, and a 
subharmonic amplitude pressure at the Hopf bifurcation 
point of |p lffB | = 0.758 MPa. From this discussion it 
follows that, in order to obtain experimentally the desta- 
bilisation of the subharmonic field, the resonator must 
be designed with a mode distribution that allows such 
values of detunings. 



Equations (|3J) have been numerically integrated by 
means of a fourth-order Runge-Kutta algorithm, and us- 
ing the acoustical parameters discussed in the previous 
section. According to this, we have chosen the numerical 
values A = 1, A x = -2.4, 7 = 7i = 523.6 (s" 1 ) and 
(3 = 0.0012 (s _1 Pa _1 ). As the control parameter we con- 
sider now the physical variable of the pump amplitude in- 
cident on the resonator \pi n \, related to the pump value 
\E\ appearing in the Eqs. © as \p in \ = (2L/c)y / l/T 
\E\. The resonator has a length of L = 10 cm, and the 
end walls a reflectivity coefficient 1Z = 0.9, which is in 
agreement with the decay rate values considered above. 

The numerical results are summarized in Fig. 4. The 
show the existence of a bifurcation of the subharmonic 
solution Eq. JBJ at the injected amplitude \pi n \ = 0.60 
MPa, as predicted by the linear stability analysis. How- 
ever, the temporal evolution of the fields near the bifur- 
cation point is not oscillatory but chaotic, as shown in 
Fig. 4(a). At the right, the corresponding phase portrait 
is shown, where the dense structure of the attractor is 
appreciated, as corresponds to chaotic behaviour. When 
decreasing the injected amplitude below the bifurcation 
point, the temporal behaviour persists, denoting that the 
Hopf instability is subcritical. As pumping is decreased 
, the temporal evolution of the amplitude follows one of 
the universal scenarios leading to chaos, namely a period- 
doubling scenario or Feigenbaum route to chaos. Since 
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P4, P2 and Pi states are observed. In the plot a grey 
scale has been used, representing pressure leves ranging 
from to 100 dB. 
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FIG. 5: Spectrum of solutions as a function of the injected 
pump amplitude, for the parameters given in Fig. 4. 



VII. CONCLUSIONS. 



FIG. 4: Time evolution (left) and phase portrait (right) of 
the sunharmonic field for A = l,Ai = -2.4, /3 = 0.0012, 
70 = 7! = 523.6, and different drivings, (a) \pi n \ ~ 0.60 
MPa. (b) |p in | = 0.58 MPa. (c) \p ln \ = 0.55 MPa. 

we are decreasing the control parameter, we observe the 
transition from chaos to periodicity in an inverted way. 
Slightly below the bifurcation point, regular oscillations 
are observed with period eight (Pg) and period four (P4), 
which exist in a narrow window. At \pi n \ — 0.58 MPa, the 
evolution is bi-periodic (P2) as shown in Fig. 4(b). For 
\pin\ = 0.55 MPa the evolution is periodic with a single 
period (Pi). This case is represented in Fig. 4(c). Fur- 
ther decreasing the control parameter, the system again 
enters in a chaotic regime, mediated a reversed period 
doubling cascade. This chaotic regime exist only in a 
narrow window, and the scenario reported above repeats 
until the value \pi n \ = 0.53 MPa, where the subharmonic 
field jumps again to the stable stationary branch given 
by Eqs. |JB}. This corresponds to the turning point of 
the dynamic branch, associated to a fold or saddle node 
bifurcation. 

The complete scenario of bifurcations can be repre- 
sented in a single plot (Fig. 5), where the frequency spec- 
trum of solutions is given as a function of the control pa- 
rameter. The subcritical character of the bifurcation, and 
the complexity of the temporal evolution in this system is 
evident here. Several chaotic windows are observed (dark 
horizontal bands, corresponding to a dense spectrum), al- 
ways followed by a period-doubling sequence where Ps, 



In this paper, the problem of subharmonic generation 
in a nonlinear acoustic resonator is treated from the point 
of view of nonlinear dynamics. A linear stability analysis 
reveals the possibility of time dependent solutions, which 
depending on parameters can be regular or chaotic. Nu- 
merical analysis are consistent with the theoretical pre- 
dictions. 

Althought the proposed model presents many analo- 
gies with similar systems studied in nonlinear optics, such 
as the two-level laser or the optical parametric oscillator, 
the peculiarities of acoustical resonators (specially the 
non-equidistant mode distribution, and the loss mecha- 
nisms) motivates a careful analysis. We conclude that 
the observability of the phenomena is fundamentally re- 
stricted by existence of low (experimentally achievable) 
threshold values, a situation which strongly depends on 
the resonator linear properties, such as mode distribution 
and quality factors. We believe that the proposed model 
could explain some experimental observations of complex 
dynamical behaviour in acoustic resonators reported in 
the past pdj . 
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